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ABSTRACT 

The structure function of opacity fluctuations is a useful statistical tool to study tiny scale structures 
of neutral hydrogen. Here we present high resolution observation of H i absorption towards 3C 138, 
and estimate the structure function of opacity fluctuations from the combined VLA, MERLIN and 
VLBA data. The angular scales probed in this work are ^ 10 — 200 milliarcsec (about 5 — 100 AU). 
The structure function in this range is found to be well represented by a power law St-{x) ~ x^ with 
index (3 ^ 0.33 ± 0.07 corresponding to a power spectrum Pt{U) ~ U~'^'^^. This is slightly shallower 
than the earlier reported power law index of ^ 2.5 — 3.0 at ~ 1000 AU to few pc scales. The amplitude 
of the derived structure function is a factor of ^ 20 — 60 times higher than the extrapolated amplitude 
from observation of Cas A at larger scales. On the other hand, extrapolating the AU scale structure 
function for 3C 138 predicts the observed structure function for Cas A at the pc scale correctly. 
These results clearly establish that the atomic gas has significantly more structures in AU scales than 
expected from earlier pc scale observations. Some plausible reasons are identified and discussed here 
to explain these results. The observational evidence of a shallower slope and the presence of rich small 
scale structures may have implications for the current understanding of the interstellar turbulence. 
Subject headings: ISM: atoms — ISM: general — ISM: structure — radio lines: ISM — turbulence 



1. INTRODUCTION 

The interstellar medium (ISM) is known to have 
clumpy density and velocity structures. Scale-free 
intensity fluctuations in a variety of tracers (e.g. H i 
21 cm emission, dust emission) is detectable over a wide 
range of scales. This scale-free fluctuation is generally 
interpreted as the signature of a turbulent ISM. The 
evidence for small scale structures and turbulence in 
the atomic ISM of the Galaxy have been found through 
various observational techniques. Power law scaling for 
the velocity dispersion of the H i 21 cm e mission has 
been observed to sub-pc scales (|Rov et al.l [20081 . On 
somewhat larger scales, the intensity fluctuations in 
H I 21 cm emission show a scale f ree p o wer spectrum 
jCrovisie r fc DickevI [TOSl iGreenI [19931 : iDickev et al] 
[2001). The H I opacity towards Cassiopeia A and 
Cygnus A show fluctuations on scales ranging from 500 
AU to 4 pc for the absorption i n the Perseus arm, the 
local arm and the Outer arm (jDeshpande et al.l 120001 : 
iRov et all 120101 ). Very long baseline interferometry 
(VLBI) obse rvations revea l struc t ures on scales as small 
as 25 AU jlpieter et al.' 1976; Diamond et al. '198'" 
Davis et al] [l996: Faison et al. 1998: Faison fc Cos; 



iStanimirovic et aD (|2010t ) have shown H i opacity 
variations on scales of 5 — 100 AU from multi- 
epoch observ at ions of some high velocity pulsars. 
iDhawan et al.l (|2000f ) have reported similar opacity 
variations on 40 — 1000 AU scales towards the mi- 
croquasar GRS 1915 -I- 105. Observations of some 
other p ulsars, however , do ii o t show any opacity vari- 
ations ([Johnston et al.l 120031: IStanimirovic et al.l [20031 : 
[Minter et al.[[2005[ : [Stanimirovic et al.|[2010D . indicating 
that such tiny scale structures may be rare. 

Small scale ISM structures are not unique to the Milky 
Way. The H i emission from the Large Magellanic Cloud, 
the Small Magellanic Cloud, several dwarf galaxies in 
the M81 group and a sample of dwarf and nearby spi- 
ral galaxies a lso show considerab le small scale intensity 
fluctuations ([Stanimirovic et al.l 1999; Westpfahl et al. 
1998 : IStanimirovic fc Lazar ian [2001 : Elmegreen et al. 
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Begum ct al. 

From an observational point of view, the H i 21 cm 
opacity fluctuations depend on fluctuations in the den- 
sity and temperature of the gas. But, for small den- 
sity fluctuations, the dep endence on temperature is weak 
([Deshpande et all l2000f l. Earlier, these small fluctua- 
tions were believed to be physical structures with den- 
sities ~ 10^ — 10^ cm~^. However, explaining the 
existence of these structures in equilibrium with or- 
ders of magnitude lower density gas in the surround- 
ing medium, is challenging. Later, it was suggested 
that the fluctuations at larger spatial scales could, in 
projection, give ri se to the observed small scale trans- 
verse fluctuations ([Deshpande 2000 ). On the other hand, 
there are some direct observational evidence of ~ 3000 
AU size H i structures with a density of ^ 100 cm~ ^ 
([Braun fc Kanekaj [20051 : IStanimirovic fc HeilesI [20051 ). 
Recent numerical simulations of the turbulent ISM have 



also indicated the existence of tiny scale structures 
(iNagashima et al] 120031: iVazauez-Semadeni et alj [20061 : 
iHennebelle fc Auditll2007l) . The short evaDoration time- 
scale (~ 1 Myr) of these structures imply that they can 
survive if either the ambient pressure is much higher than 
the average ISM pressure or they are formed on a similar 
time-scale. In any case, their physical properties are not 
yet well understood. 

Here we present results on tiny scale H i opacity fluc- 
tuations towards the quasar 3C 138 derived from VLBA, 
MERLIN and VLA data. Combining MERLIN and VLA 
data with VLBA data has resulted in improved sensitiv- 
ity and larger spatial dynamic range to detect small op- 
tical depth fluctuations and to derive the structure func- 
tion. The data analysis method is outlined in Section Sj2l 
Section [J3] contains the results, and the conclusions are 
presented in Section [|4l 

2. DATA ANALYSIS 

For this analysis, we have used VLBA, MERLIN and 
VLA data. Results from the MERLIN and VLBA 
observ a tions were pr eviously presented by Davis et al. 
(jHH), iFaison et all (jB98) and Brogan ct al. (2005). 
The VLBA observation (Project BD0026) was carried 
out on September 10, 1995. The spectral resolution of 
this data is -^ 0.4 km s"^ The MERLIN observation 
was carried out on October 22 and 23, 1993, and the 
spectral resolution in this case is ~ 0.4 km s~^. The 
VLA A-configuration observations were carried out on 
September 07 and 13, 1991 (Project AS0410 and TEST). 
The channel width and spectral resolutio n in this case 
are ~ 0.8 and ~ 0.6 km s~^, respectively. iBrogan et al.l 
()2005( ) reported temporal variation of H i opacity towards 
3C 138 from the VLBA observations taken in 1995, but 
also include data from 1999 and 2002. Here, we have only 
included the 1995 data, since it is closest to the MER- 
LIN and VLA observation epochs. The longest base- 
line was about 0.18, 1.0 and 36 MA (angular resolution 
of about 1000, 200 and 6 mas) for the VLA, MERLIN 
and VLBA data respectively. However, even with the 
combination of the VLBA, MERLIN and VLA, sufficient 
brightness temperature sensitivity was achieved only on 
angular scales larger than about 20 mas, corresponding 
to baselines shorter than about ~ 15 MA. 

Data analysis was carried out using the NRAO As- 
tronomical Image Processing System (AIPS). All the 
data were converted to J2000.0 epoch, and the spec- 
tral axis was regrided to a common spectral resolution 
of 0.4 km s~^. A combined dataset was produced af- 
ter first reweighting the VLBA, MERLIN and VLA data 
to have similar data weights (i.e. so they would con- 
tribute equally at the imaging stage). A continuum and 
continuum-free dataset were produced by fitting the line- 
free channels in the uv-plane. Self-calibration was re- 
peatedly done on the continuum data starting with a 
VLBA image as the initial model, and then including 
the MERLIN and the VLA baselines at later stages. The 
self-calibration solutions obtained for the continuum data 
were also applied to the spectral line data cube. The final 
continuum image, shown in Figure ([T|), and the spectral 
line image cube were made using multiscale imaging with 
a 20 mas restoring beam. The structures of the contin- 
uum and line images agree with the VLBA images of 
IBrogan et al.l ()2005[ ). The spectral image cube and the 
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Figure 1. 1.4 GHz continuum image of 3C 138 from the combined 
VLBA, MERLIN and VLA data with a 20 mas restoring beam. 
The rms noise of the image is 12 mjy/beam. The contour levels 
at (1.0, 1.4, 2.0, 2.8, 4.0, 5.6, 8.0, 11.2, 16.0, 16.8) x36 mJy/beam 
are overplottcd on the gray scale image with linear intensity scale. 



continuum image were used to obtain the optical depth 
image cube. A 3ct cutoff was used for blanking the contin- 
uum image, and the same blanking window was applied 
to the spectral data to get rid of possible contribution 
from high noise or other low level image artifacts. The 
optical depth image cube was used to derive the opacity 
fluctuation structure function. 

We note that the structure function amplitude of the 
line-free channels is more than an order of magnitude 
lower than the structure function amplitude of the chan- 
nels with H I absorption signal. But, for the channels 
with H I signal, the optical depth noise is significantly 
higher than that of the line-free channels. Since the rms 
noise in the continuum image (trp) is significantly lower 
than the rms noise in line channels {(Tj), the optical depth 
noise can be approximated as cr,- « e'^{ai/Io), where Iq 
is the continuum image flux density. Thus, for high opti- 
cal depth, the noise is scaled up by a factor of e"^. So, the 
optical depth noise is highly correlated with the optical 
depth itself. Numerical analysis shows that, as a result, 
the noise bias will have a similar power law dependence 
on all angular scales, and will scale up the structure func- 
tion amplitude without changing the power law index. 
Moreover, the H i brightness temperature for the line 
channels will contribute to the system temperature, and 
will also increase the image noise. As a result, estimat- 
ing the noise contribution from the line- free channels will 
underestimate the errors to the structure function due to 
optical depth uncertainties. 

Statistical errors on the derived structure function for 
individual channels are estimated from the variance of 
the measurements within the bin, and do not account for 
the high optical depth measurement uncertainty. If the 
high opacity pixels are excluded, the effect of this uncer- 
tainty and the noise bias to the derived structure function 
is minimized. But in that case, the opacity image is mod- 
ified by a window function, and the structure function 
can only be derived for about a factor of 4 smaller range 
of angular scales, and only ~ 8 — 10 spectral channels. 
Hence, first the high opacity pixels are included to derive 
the initial structure function. With significant noise bias, 
the amplitude of the structure function in this case is only 
an upper limit. However, one can use a larger range of 
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Figure 2. Integrated H I 21 cm absorption spectrum towards 
3C 138 from combined VLB A, MERLIN and VLA data. The chan- 
nels in the shaded region are used to derive the structure function. 
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Figure 3. Structure function from a few selected channels with H I 
absorption (dashed line) , and from one channel with no absorption 
(solid line). For clarity, error bars are shown for only one channel. 



angular scale to fit the power law index. Then t j Or cut- 
off is used to minimize the noise bias and to compute the 
amplitude of the structure function from a limited range 
of angular scales. For both cases, statistical errors are 
estimated using the channel to channel variation of the 
structure function. Essentially, the opacity values in any 
channel may be considered as an independent realization 
drawn from a parent distribution of opacities with asso- 
ciated uncertainty. So, the channel to channel variation 
of the structure function will include the opacity uncer- 
tainty effect. The errors are estimated assuming that 
there is no spectral variation of the structure function. 
Thus, one can only derive an average structure function 
from all the spectral channels. 

3. RESULTS 

3C 138 is unresolved at the short baselines (< 45kA) 
of the VLA. Thus, by combining the VLBA, MERLIN 
and VLA data, we are able to image 3C 138 maintain- 
ing the small scale structures at the VLBA resolution, 
and, at the same time, restoring the total flux density 
observed by the VLA. The intermediate MERLIN and 
VLBA baselines provide a wide range of uv-coverage, 
allowing us to probe structures at scales ranging from 
about ten milliarcsec to about 0.20 arcse c. Adopting 
a dis tance to the absorbing gas of 500 pc (jFaison et al.l 
Il998f ). the range of angular scales translates to 5 — 100 
AU. The integrated H i 21 cm absorption spectrum is 
shown in Figure ([2]). The shaded region in Figure ([2]) in- 
dicates the channels that have adequate signal to noise in 
optical depth images for further analysis of small scale 
opacity fluctuations. We have calculated the structure 
function for 36 individual spectral channels, each with a 
velocity width of ^ 0.4 km s"""^, and used all the mea- 
surements together to derive the best estimate of the 
structure function slope and amplitude. 

For quantitative study of small scale fluctuations, sta- 
tistical tools like the correlation function, structure func- 
tion, and power spectrum are very useful. The observed 
shape of these functions are expected to be related to 
the physical processes that gives rise to the fluctuations. 
The power spectrum of the opacity fluctuations PriU) is 



the Fourier transform of the correlation function ^t{x): 
Pr{U) = Pr{u,v) ^ f /er(^m)e-2"("'+'"")dWm, 

J J (1) 

where (l, m) and (u, v) are the sky coordinate and inverse 
angular separation respectively. Here, ^t{x) is defined as 

^7-(a;) = ^tC — ^',"T- — Ti') = ((5r(/, m)(5T(/', m')) , (2) 

and the structure function Sr{x) — [At(x)]^ is given by 

Sr{x) =Sr{l~ I', m - m') = ([t(/, m) - t{1\ m')f) . (3) 

Strictly speaking, Pt and ^7. are function of U and x, 
respectively, but we have assumed statistical isotropy 
for the purpose of simplification. Thus, throughout this 
analysis, Pt and ^t- are considered to only be func- 
tions of scalar amplitude U and x, respectively. Gen- 
erally, if the power spectrum is a power law of the form 
Pt{U) = PoU~°' with 2 < a < 4, the structure function 
St{x) will have a power law index of a — 2, and the rms 
opacity fiuctua tion AtJx) will have a power l aw index of 
(a ~ 2')/2 (Lee fc Jokippiil [TqtI IDeshpandel [200Q1 . We 
also note that if the density fluctuations are small (that 
is Sp/{p) << 1), the density fluctuation power spectrum 
is expected to be very similar to the opa city fluctuation 
power spectrum (jDeshpande et al.ll2000( ). 

The derived structure function of the H i absorption 
towards 3C 138 for a few spectral channels are shown in 
Figure ([3]) . The amplitude and the shape of the function 
changes from channel to channel. The structure function 
from one spectral channel with no H i signal is also shown 
in this figure for a comparison. The measured struc- 
ture function for all 36 spectral channels with adequate 
signal to noise is shown in Figure (J4l). The open point 
symbols are structure function values without any t/ctt- 
cutoff. Clearly, below ~ 10 milliarcsec, the structures 
are smoothed by the synthesized beam. Above ~ 200 
milliarcsec, the effects of the spatial window due to the 
shape of the continuum emission significantly modify and 
flatten the measured structure function. The mean and 
rms values of these points in each bin are shown as filled 
circles with error bars (with position shifted slightly to 
the left for clarity). The best fit power law for the ensem- 
ble (dash-dot lines in the figure) has a power law index of 



0.33 ±0.07 (3(t), which corresponds to a power spectrum 
with a power law index a — 2.33. 

Next, we have used a cutoff based on t/ctt- to minimize 
the noise bias and obtain a better estimate of the struc- 
ture function amphtude. As mentioned earher, with this 
cutoff we have a reUable estimate of the structure func- 
tion only for ^ 8 — 10 channels and a limited range of 
angular scales. Here, we have used a cutoff at t/cti- = 2 
with usable range of ~ 10 — 40 mas for eight spectral 
channels. These estimates (shown as filled points joined 
by dotted lines in the same figure) are also consistent 
with a power law index of 0.33. Since r/cr,- > 2, the sig- 
nal structure function is at least a factor of 4 higher than 
the noise bias case. So, the amplitude of the signal struc- 
ture function is constrained to be higher than 0.8 times 
the observed amplitude. Unfortunately, for any higher 
cutoff, the flattening due to the spatial window domi- 
nates, and essentially makes it impossible to extract the 
true structure function. The solid line in th e figure is the 
structure function from iDeshpandd ()2000D extrapolated 
from the observed structure function power index of 0.75 
in the range of ~ 0.02 — 4.0 pc. Clearly the observed 
structure function amplitude at smaller scales (5 — 100 
AU) is ~ 20 — 60 times higher than the extrapolated 
value. In terms of rms opacity fluctuation, the observed 
values are a factor of 4— 8 times higher than the predicted 
values. 

Apart from considering measurements from all the 
channels together for fitting, each channel is considered 
separately to cross-check if there is any significant chan- 
nel to channel variation of the structure function. Note 
that the errors on the structure function derived from sin- 
gle channel refiect only the variance of the measurements, 
and do not include the effect of large optical depth uncer- 
tainties. The best fit power law amplitude and index for 
each channel are shown in Figure ([SJ . Error bars for indi- 
vidual channels are 3cr error from the power law fit. The 
structure function amplitude is a factor of 2 — 5 higher 
for channels with stronger absorption. This is consistent 
with the fact that the strong noise bias in these chan- 
nels are expected to scale up the amplitude. However, 
we note that the average value of the power law index 
0.36 ± 0.06 (3cr interval shown in dotted lines) is consis- 
tent with no channel to channel variation of the index. 
This best fit value of the power law index is also consis- 
tent (at better than 3cr uncertainty level) with the power 
law index value of 0.33 ± 0.07 we derived earlier from the 
ensemble of structure functions from all the channels. 

Typically, the power spectra from the H i emis- 
sion and absorption observations have a power law 

index a ~ 2.5 — 3.0 at scales of ^ 1000 

AU t o a few pc (ICrovisier fc Dickcvl 119831: IGreenI 
1993'; Stanimirovic et al. 1999; Dcshpandc et al. 2000: 
Begum et al.. .2006. : .Dutta et al.. ■2009a .b.: .Rov et al.. 
20101 ). For 3C 138, the derived power law index (a ~ 
2.33) for the power spectrum is slightly shallower. This 
flattening may arise from two reasons. Earlier it has been 
reported that, in the case of turbulence in a supernova 
remnant, a transition from 3-d to 2-d turbulen ce makes 
the power spectrum shallower (jRov et al.ll2009t ). So, in- 
stead of homogeneous random fluctuations, anisotropic 
turbulence at smaller scales can result in a shallower 
spectrum. These fluctuations may be averaged out for 
large scale observations with an arcsec beam. The other 



plausible reason may be small scale physical processes 
like viscous damping. Viscous damping for incompress- 
ible magnetized turbulence can make the mag netic en- 
ergy spectrum significantly less steep (|Cho et al. 2002). 
Since, even with very small ionization fraction, H i is 
coupled to the Galactic magnetic field, any possible cor- 
relation of H I density with the magnetic field strength 
will also flatten the opacity fluctuation power spectrum. 
In this regard, it will be interesting to check if simula- 
tion of compressible magnetohydrodynamic turbulence 
including damping and anisotropy at small scales can 
explain these observational results quantitatively. 

Our results also point out an amplitude mismatch be- 
tween the observed structure function and that of the 
structure function extrapolated from large scale. The 
derived amplitude is more than an order of magnitude 
higher at tens of AU scale for 3C 138. The observed 
higher amplitude at small scales may be due to few pos- 
sible reasons. First of all, for our analysis, we have used 
3C 138 (I = 187°, b = —11°) a s a background source, 
whereas Deshpande et all (|2000[ ) used Gas A (/ = 112°, 
b = —2°) for deriving the structure function at pc scales. 
These two lines of sight are widely separated, probing 
H I in different parts of the Galaxy, where the small 
scale structures may, in principle, be intrinsically differ- 
ent. Alternatively, it may be an effect of the uncertainty 
in the measurement of structure function towards Gas A, 
and extrapolation over a large angular scale giving rise 
to this mismatch. In this regard, we find that extrapo- 
lating the shallower structure function for 3C 138 to 0.5 
pc scale predicts the large scale structure function am- 
plitude within ~ 10% of the observed value for Gas A. In 
other words, the amplitude of the structure function at 
AU scales is, as expected in this second possibility, com- 
pletely consistent with the observed rms opacity fiuctua- 
tions at the pc scale. Finally, it is possible that there are 
physical mechanism(s) which modify both the amplitude 
and slope of the structure function at AU scales. For 
example density fluctuations due to intermittency may 
in part ex plain this excess of tin y scale structures. We 
note that, (Fa lgarone et al.l l 20 09') have reported millipar- 
sec scale extreme velocity shears, which are explained as 
signature of intermittency. Further observations for more 
lines of sight, and measurements of structure function to 
bridge the gap between AU and pc scales will be very 
helpful in distinguishing between these different possible 
scenarios. 

4. GONGLUSIONS 

We have used high resolution images of H i absorp- 
tion towards 3C 138 from the combined VLA, MER- 
LIN and VLBA data to study the small scale opacity 
fluctuations. The fluctuations have a power law struc- 
ture function with an index of /3 ^ 0.33 ± 0.07 over the 
scales of '--^ 5 — 100 AU. This implies a power spectrum 
Pr{U) ^ U^^'^^, slightly shallower than the earlier re- 
ported values of a ~ 2.5 — 3.0 from observations over 
larger scales. The amplitude of the structure function at 
the AU scale is found to be significantly higher than the 
extrapolated amplitude from observations of Gas A at pc 
scale. The observed rms opacity difference for 3G 138 is 
about 4 — 8 times higher than the prediction from this 
extrapolation from Gas A. However, extrapolating the 
observed AU scale structure function for 3G 138 consis- 
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Figure 4. Measured opacity fluctuation structure function from all spectral channels with adequate signal to noise. Open squares are 
estimates without a r/crr cutoff. The mean and rms values in each bin are shown as filled circles with error bars (with position shifted 
slightly to the left for clarity). The best fit power law (shown as dash-dot lines) for the ensemble has a power law index of 0.33 ± 0.07 {3a). 
The structure function derived using t/cTt- = 2 cutoff (filled squares joined by d otted lines) is also consistent with a power law index of 
0.33. The solid line is the extrapolated structure function from Deshpandcj I I2000I) . 
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a shallower slope and the presence of rich structures at 
smaller scales. Further measurements of structure func- 
tion between AU and pc scales, and detailed numerical 
simulations to quantitatively verify these results can im- 
prove the understanding of the interstellar turbulence. 
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